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In a number of experimental situations, single polymer molecules can be suspended in a vacuum. 
Here collisions between such molecules are considered. The limit of high collision velocity is inves- 
tigated numerically for a variety of conditions. The distribution of contact times, scattering angles, 
and final velocities are analyzed. In this limit, self avoiding chains are found to become highly 
stretched as they collide with each other, and have a distribution of scattering times that depends 
on the scattering angle. The velocity of the molecules after the collisions is similar to predictions of 
a model assuming thermal equilibration of molecules during the collision. The most important dif- 
ference is a significant subset of molecules that inelastically scatter but do not substantially change 
direction. 

PACS numbers: 



I. INTRODUCTION 

Although polymer molecules are most commonly stud- 
ied in solution or in solid form [l[ , there has been increas- 
ing technological use for them in a vacuum. They are 
often prepared in this state as part of the technique used 
to identify protein molecules with mass spectrometry Q . 

Recently the properties of single chain molecules in a 
vacuum were studied theoretically and by means of com- 
puter simulation [3-0| . It was shown that such molecules 
have unusual statistical properties and that the dynamics 
are very different from those found for similar molecules 
in solution. With no excluded volume, the lack of a sol- 
vent means that the only damping that can occur is in- 
ternal to the chain and this leads to slowly damped os- 
cillatory behavior for time dependent correlation func- 
tions of polymer position. With excluded volume, the 
time constant for relaxation scales with chain length N 
as 7V^-^^^ °^ substantially faster than for corresponding 
chains in solution [l|. Oscillatory behavior is quite pro- 
nounced for short chains f6| but is suppressed when they 
are longer 3] . 

The equilibrium size of such chains is also influenced by 
the conservation of angular momentum Q , and the exact 
statistical properties of an ideal chain with this constraint 
can be calculated. When the total angular momentum is 
zero, the radius of gyration is smaller relative to a chain 
without this constraint enforced. 

In experimental situations, these polymer chains are 
often charged and are accelerated by an external field. 
As far as the author is aware, there have been no exper- 
iments that have measured their conformations in this 
state. Such experiments would be very interesting and 
might allow for the probing of additional features that 
could be used to characterize the chemical composition 
of such chains. An important aspect in understanding 
this situation are the nature of collisions between chains. 

In equilibrium, collisions will obey detailed balance and 
the average properties of two chains before and after a 
collision will be identical. Under the conditions neces- 
sary to perform mass spectrometry, nonequilibrium con- 



siderations become necessary. Because the molecules are 
charged and accelerated by large fields, they acquire a 
high center of mass velocity. It seems likely that under 
such circumstances, one would have collisions occurring 
between molecules with high relative velocities, as their 
charges and masses are not all identical. Such collisions 
differ from most collisions in that they are occurring be- 
tween very large molecules, and are expected to have 
quite a different character than previously studied. 

In addition it should be possible to directly study colli- 
sions between molecules with high relative center of mass 
velocities using modifications to the highly sophisticated 
apparatus that is in current use such as matrix-assisted 
laser desorption/ionization mass spectrometry (MALDI 
MS). The acceleration voltage typically of order order 
lOkV 8]. So with a single charge on a protein this 
amounts to a center of mass kinetic energy of W^eV. 
Compared to the center of mass thermal energy at AOOK, 
this over 10^ times greater. It should be noted that 
this large field does not cause further ionization of the 
molecules because the acceleration occurs over centime- 
ters and so the electric field is still small compared with 
that needed for ionization. There are many sophisticated 
variations [9] of MALDI that suggest that modification 
of the apparatus for the purpose of studying collisions 
should be feasible. 

Therefore it is of interest to explore situations where 
the collision velocities of molecules are much larger 
than the thermal values given by the equipartition theo- 
rem [13] and this is done in this paper. 

Throughout this work we will consider collisions in the 
center of mass frame. Transforming to other frames is 
straightforward. For highly inelastic collisions, such as 
the ones described below, this is the most natural refer- 
ence frame to use. 



II. THE HIGH COLLISION VELOCITY LIMIT 

Self avoiding chains are characterized by a mass per 
monomer m, step length /, and an excluded volume pa- 
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rameter. Excluded volume can be understood as an ef- 
fective hard core radius which prohibits monomers from 
getting closer than a certain distance. For the sake of 
generality, consider that the chains have different num- 
bers of monomers Ni and and corresponding masses 
Ml — mNi and M2 — mN2. We will consider the col- 
lision in the center of mass reference frame, so that ini- 
tially, the total momenta of the molecules are P/ and 
— P/. In the limit where Pj is very large and the in- 
ternal energy of a chain is kept constant, the internal 
energy of a chain becomes negligible in comparison with 
the center of mass kinetic energy. During a collision, we 
will see that typically, a substantial fraction of this mo- 
mentum is transfered into internal kinetic energy of the 
chains. Therefore in this limit, the internal energy of a 
chain before a collision can be neglected, and we can set 
the initial kinetic energy of the chains equal to zero. In 
this case, with only hard core potentials, the scattering 
of self avoiding chains only depends on velocity as a pref- 
actor. The angular dependence of scattering becomes 
independent of P/, and the distribution of the final scat- 
tered momentum Ppi only depends on Pp/Pj. This is 
an interesting limit to consider because of the lack of de- 
pendence on the chains' temperature, and therefore this 
will be studied in detail below. 

The distribution of resulting directions emerging from 
the collision is an interesting quantity to examine. The 
angle of deflection 9 (the inclination angle) after a col- 
lision can be anisotropic, but the azimuthal angle 
must always be isotropic. We use the convention that 
7r/2>6'>0is the direction of forward scattering. This 
is shown in Fig. [1] 
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A. Thermal Limit 

When the chains are in contact for long enough that 
they have fully equilibrated, the energy in the velocity 
degrees of freedom will then be described by the equipar- 
tition theorem with the caveat that they strictly obey 
conservation of energy, momentum, and angular momen- 
tum. For a large number of monomers, momentum and 
angular momentum conservation make a negligible cor- 
rection to the energy in each degree of freedom. The 
initial energy is 



E=—Pi^{— + —) 



(1) 



Assuming thermal equilibrium, when the two chains col- 
lide, we can define a temperature T, and statistical me- 
chanics gives the relationship between this and the en- 
ergy E. The average total kinetic energy is {K) — 
[Ni -f N2)dfkBT/2, where df are the number of degrees 
of freedom per monomer, so in the special case of an 
athermal system, E = K so that 
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We would like to calculate the average center of mass 
energy of both molecules after the collision. The posi- 
tion and momentum degrees of freedom are independent 
so we need only consider kinetic energy corresponding to 
particles with momenta Pj,i, and Pj^2 where the first sub- 
script, j, indexes the monomer, and the second subscript 
indexes the polymer. 
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The probability of finding the system with a given set of 
momenta is proportional to exp{—/3K). subject to the 
constraint the 
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FIG. 1. (Color Online) A collision in the center of mass 
frame. Polymers A and B move toward each other as shown 
by the solid line. After the collision they come apart as shown 
by the dashed lines. The angle 6 is the angle of deflection. 



When two chains collide, they remain in contact for 
some period of time, tc, that depends on the precise de- 
tails of the initial configurations of the chains. Because 
the chains are flexible, we expect to have highly inelas- 
tic collisions. In the limit of a very long tc, much longer 
than the relaxation of a chain, the two chains will have 
fully thermalized. We will now describe what we expect 
in this limit. 



We would like to find the probability distribution of the 
total momemtum of one chain 



Pi 



Ni 



(5) 



This problem is mathematically identical to the problem 
of a gaussian ring polymer with Ni + N2 monomers. In 
particular, one can regard the particle momenta as dis- 
placement vectors between nearest neighbors. Then Pi 
can be interpreted as the displacement of two monomers 
separated by A^i monomers. For the purposes of this cal- 
culation, this is the same as two linear chains in parallel. 
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of lengths Ni and N2 . Therefore probability distribution 
of Pi is 



Therefore 



P(Pi) oc e ^ U./i + M2^ 



(6) 



(7) 



where l/M^ — 1/Mi + I/M2 is the reduced mass. 

Alternatively this result can be seen by considering 
a gas in thermal equilibrium and considering the joint 
momentum distribution of two molecules of mass Mi and 
M2, P(pi,p2) oc exp{—f3K). We are interested in the 
distribution of momentum in the center of mass reference 
frame. The kinetic energy K is the sum of the internal 
plus center of mass kinetic energy. Because the internal 
momenta are equal and opposite, the internal term Ki is 



2 ^Afi M2 



(8) 



and is independent of the center of mass momentum. 
This implies the distribution of Pi is of the form of Eq. 

El 

In the athermal limit, we can use the value of T given 
by Eq. [5] and substitute that into Eq. [7] obtaining 
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(9) 



Therefore the variance of the final center of mass speed 
of chain 1, Vcm,i,F, is related to its initial speed Vcm,i,i 
as 
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cm.l,/ 



df{Ni + N2) 



(10) 



In a gas of molecules with a large number of internal 
degrees of freedom, the distribution of total molecular 
energy is highly peaked. The standard deviation of this 
energy is smaller than the mean by a factor of 1/ \/N. 
In the limit of large N we can ignore ignore fiuctuations 
in this energy. In thermal equilibrium, we would like to 
know the distribution of center of mass speeds seen that 
result from a collision. 

To derive this distribution, assume we have just two 
molecules in a container with periodic boundary condi- 
tions with total center of mass momentum of zero. As- 
sume the size of the container is much larger than the 
size of a molecule and that the system is in thermal equi- 
librium. Such a system is expected to be ergodic and 
hence be able to reach thermal equilibrium. The dis- 
tribution of the center of mass speed measured at one 
time is Maxwellian. This differs from the distribution of 
speeds that are measured after a collision. In the latter 
case, we are not measuring speeds at arbitrary times, but 
rather, after a collision. This difference in sampling al- 
ters the distribution. After a collision the center of mass 



of molecule 1, will move in a straight line with speed 
Vcm,i.F until it suffers another collision. The time spent 
in this state before undergoing another collision will be 
oc l/fcm,i,F- So to compenstate for this, the distribution 
of final speeds must be multiplied by an extra factor of 
Vcm.i F so that 



P{Vcm,l,F)dVcm,l,F OC vl^ i pt 
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Here cr^ is the variance of one of the components of the 



final velocity, and therefore cr^ 
Eq. [H 
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For a gas in thermal equilibrium, the distribution of 
speeds of molecules striking a wall is well known from 
the study of effusion [ll| and this is the form that is 
expected in that case as well. 

The distribution of directions for the final velocity in 
this limit should be isotropic. The only subtlety being 
that the total angular momentum Lt is conserved, so if Lt 
is zero, it might be thought that this would preclude the 
two chains rotating relative to each other and that this 
would cause them separate in a way that preserves their 
initial relative orientation. However it is still possible to 
rotate the two chains by circular motion of only parts 
of them. For example, if an end monomer rotates by 
one revolution, this will cause the rest of the system to 
compensate by rotating in the opposite direction. This 
will cause a rotation in the relative positions of the two 
chains, without changing the total angular momentum. 
Therefore in the limit of long collision times, the relative 
orientation of the two chains is not conserved and the 
final direction should be isotropic. 

In terms of the variable 



= cos 



(13) 



the distribution of z should be uniform. For this reason 
we will use z rather than 9 to characterize the results of 
the simulations discussed below. 

Note that the above consideration could apply to any 
large molecules undergoing collisions, not just polymers. 
The only thing required is that they remain in contact for 
sufficiently long that thermal equilibrium is established. 



B. Simulation 

We model self avoiding chains as was done previ- 
ously [3]. The distance between links is maintained at 
a constant value of Z = 1. The chain is modeled as freely 
hinged and there is no chain stiffness. A repulsive poten- 
tial between all monomers was included. If the distance 
between two monomers is r, their potential was taken to 
heVir) = 40(l-(r//)2)5. For r > I, V{r) =0. A diverg- 
ing hard core was not used for the sake of efficiency. The 
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potential at the center is very high making chain crossing 
exceedingly unlikely and no chain crossing was observed. 

Rigid links were chosen to avoid equilibration problems 
that can occur due to the quasi-one dimensional nature 
of this system [l^. The simulation method solves 
Newton's laws for this system so that all conservation 
laws are well satisfied. It also implements the rigid link 
constraints in an efficient way, using 0{N) operations for 
every integration step. The method is described in detail 
inRef. 0. 

Many collisions were studied with different initial con- 
figurations of the chains and their statistics were ana- 
lyzed. Wc took the initial center speed of chains to be 1, 
and TO = I. The time step for the simulations was 0.02. 

The chains were initially equilibrated for many relax- 
ation times. For example for N = 128 the equilibration 
time was 40, 000 steps. Then their velocities were made 
negligibly small by applying a large damping for many 
damping times so that their final velocities were less than 
10^^. These chains were given their initial center of mass 
speeds Wcm,i,/ = 1 in the direction between the two center 
of masses. This way there was no initial internal energy 
of the chains, corresponding to the limit of high relative 
collision velocity. The chains were initially separated by 
at least half of their total chain length. This ensures that 
they are well separated before a collision takes place. To 
obtain adequate statistics, many collisions, about 10'*, 
were used with different random initial conditions. 

The length of time they remained in contact was de- 
termined by monitoring how the center of mass veloc- 
ity changed. Before a collision, as the chains ap- 
proach each other, Vcm remains constant, but then as 
soon as the chains collide, there will be a change in this 
velocity. The point where Vcm first changes signals the 
start of a collision. In order for a collision to be consid- 
ered over, Vcm had to remain constant for all subsequent 
times. Time evolution was not stopped until the chains 
had completely separated from each. This was imple- 
mented by requiring that the evolution continue if there 
were two monomers from different chains that were closer 
than half their original separation at the start of the sim- 
ulation. This eliminated the cases where the two chains 
appeared to have separated but later collided again. The 
final speeds were also recorded, as was the angle of de- 
flection 9 of the chains after the collision (again, with the 
convention that 7r/2>0>Ois forward scattering). 



C. Results 

Fig. [2] shows the distribution of scattering angles as a 
function of z, as defined in Eq. [T31 For A'^ = 32, there 
is more scattering backwards, that is, in the direction 
of negative z. For larger TV, this changes and for N = 
128, the effect of backwards scattering is much weaker. 
Instead, there is a strong peak for z — 1. This means that 
a substantial fraction of collisions defiect the polymers by 
a very small amount. We will discuss this further below. 
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FIG. 2. (Color Online) The distribution of scattering angles 
9. The horizontal axis is z = cos 6. The highest curve on 
the left is for A'' — 32, the middle A'^ = 64, and the lowest 
N = 128. 



Assuming an athermal model, the temperature given 
by Eq. [5] is fc^T — 1/2. The distribution of colhsion 
times is shown in Fig. [3] for A^ = 128 (lower graph) 
and A^ = 64 (upper graph.) The average time computed 
from this distribution, for N = 128, is 210. This is com- 
parable but larger than the equilibrium relaxation time 
trei seen for similar simulations of equilibrium systems. 
The latter was determined from measuring the time de- 
pendent autocorrelation function as has been described 
previously, see Ref. [1], Fig. 3. This suggests that dur- 
ing the majority of these collisions, the system should be 
close to thermal equilibrium. There are two sets of data 
shown on each plot in Fig. [31 One where all scattering 
angles are included (triangles), the second only includes 
forward scattering, z > (squares). With only forward 
scattering, there is a prominent dip in the distribution 
at t « 150 (A^ — 128.) For longer times, the two curves 
converge quite closely. The reason for this unusual small 
t behavior will also be discussed below. 

The distribution of final velocities is shown in Fig. 
Ufa). The smooth curve without error bars, is a fit to 
the form expected in the thermal limit, Eq. [TTJ Note 
that there is a much longer tail in the simulation data. If 
instead one excludes strong forward scattering, z > 0.9, 
as shown in Fig. Hl^b), the fit to the thermal distribution 
is much improved. The fit gives = 0.064, where as the 
predicted value from Eq. [T^ the thermal limit, is 0.044. 
The distribution in Fig. IH^b) is still slightly too wide 
and this suggests that low final velocities are suppressed 
compared to the thermal prediction. 

The distribution of the maximum end to end distances 
R, occurring during scattering is shown in Fig. [5] for 
A^ = 64 and A^ = 128. There is typically substantial 
stretching that occurs during a collision. For compari- 
son, statistics for the end to end distance Re for A^ = 128 
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FIG. 3. (Color Online) The distribution of collision times t 
for a chain with TV = 128 (lower graph) and A'^ = 64 (upper 
graph), and the same distribution excluding those where the 
scattering angle is greater than tt/2. The latter distribution 
has a sizable dip at t ~ 150 for A'^ = 128. 



were calculated in thermal equilibrium for a total angu- 
lar momentum of zero. The same chain has an rms end 
to end separation of ^ (R^) = 17.5, However the max- 
imum in Fig. [5] is at i? = 38 and has a significant tail 
stretching to past R = 90. After the chains have sep- 
arated, the stretching will disappear and tend towards 
their equilibrium values. 

To summarize what we have learned so far, collisions 
between self avoiding chains in a vacuum appear be quite 
close to what one would expect in the thermal limit. For 
the longest chain studied, N = 128, the scattering is 
close to isotropic with the exception of a spike for strong 
forward scattering. The time that the chains remain in 
contact is typically longer than a relaxation time. The 
distribution of final speeds appears quite comparable to 
what one would expect from a thermal distribution ex- 
cept for collisions coming from strong forward scattering. 
Finally, at some point while the chains are in a collision, 
they are typically quite stretched. 

From this it appears that these chain collisions are 
highly inelastic and they have properties close to that of 
collisions in the thermal limit. However there are some 
interesting deviations from this limit as was noted above. 
There appears to be a fraction of collision where the angle 
of scattering is very small. If these chains are excluded, 
then the collisions have statistics much closer to that of 
the thermal limit. A simple explanation for this would be 
that this sub-population of chains are not actually collid- 



ing. However their final velocities are substantially dif- 
ferent from their initial ones, so this effect is clearly more 
subtle. By observing such strong forward scattering col- 
lisions it appears that their behavior can be understood 
as follows. 

The dynamics of chains during collisions can be under- 
stood in terms of two factors, entanglement effects and 
transfer of momentum. During strong forward scattering 
collisions, it appears from visualizations, that the chains 
do strongly collide but only weakly entangle. Different 
parts of the chains come into contact and when they do, 
these parts collide very inelastically. That is, they do not 
bounce off one another but remain in contact for some 
length of time. Therefore momentum transfer between 
the chains is still inelastic. If two inelastic masses col- 
lide, even with unequal mass, their final velocity will still 
be along the same line as the initial velocities. There- 
fore all such inelastic collisions do not change the final 
direction. Because the chains have not entangled, they 
move past each other, but have not remained in contact 
long enough to thermalize. Their final directions are al- 
most unchanged but their speeds have been substantially 
reduced due to momentum transfer. 

The above observations can be seen by examining two 
dimensional distributions, for N = 128 in Fig. [B] In Fig. 
[6l[a), the distribution of collisions is binned in terms of 
final collision speed v, and z. One can see that for strong 
forward scattering, z « 1, there are a large number of 
collisions with final speeds much larger than the thermal 
velocity. In In Fig. El^b), the distribution of collisions 
is binned as a function of collision time t and z. Again 
for strong forward scattering, there is a high peak for 
atypically short times, indicating that the collisions are 
weaker than for other scattering angles. Therefore these 
strong forward scattering collisions are short and do not 
slow down the chains all the way to their thermal values. 

The plots in Figs. [3] and S] are obtainable from these 
two dimensional distributions. The enhancement in short 
time collision frequency seen in Fig. [3] for z > is there- 
fore caused by the subset of strong forward scattering 
collisions, discussed above, that only collide for short 
times. All other collisions that more completely ther- 
malize, usually remain in contact for longer times. The 
fact that there is a dip in the distribution of Fig. [3] when 
considering collisions with z > suggests that the quasi- 
thermal chains having z > have a peak in their collision 
time distribution at i w 200. 



III. DISCUSSION 

It is of interest to examine the scaling of collisions with 
N. The radius of gyration of a swollen chain R ^ N'^ 
with « 3/5 in three dimensions. According to our as- 
sumption that the initial kinetic energy of the chains are 
zero, the probability that a monomer from chain A will 
collide with chain B is obtained by looking at the pro- 
jection of the density of chain B into a two dimensional 
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FIG. 4. (Color Online) (a) The distribution of final collision speeds v including all scattering angles and A'^ — 128. (b) The 
same distribution ignoring those where cos 6 < 0.9. For comparison the distribution for L = 64 is also shown. 
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FIG. 5. (Color Online) The distribution of the maximum end 
to end distances R, occurring during scattering for chains with 
iV = 64 and 128. 



plane. The average density of this projection is propor- 
tional to the probability of a collision. The projected den- 
sity is proportional to iV/i?^ ~ AT-i/s, There- 
fore as ^ cx) the probability that a given monomer in 
chain A will collide with chain B goes to zero. On the 
other hand, the total number of contact points in chain 
A is N^/^. Therefore for large N the molecules will suf- 
fer a collision but the above argument suggests that the 
number of monomers actually making contact will be a 
small fraction of the total. 

If the amount of time that these points remain in con- 
tact is not long enough, the chains will pass each other at 
some reduced center of mass velocities. If the collisions 
of contact points are highly inelastic, only the speed of 
the chains, but not their direction, will change. Entan- 
glements will change the above arguments considerably. 
Only one entanglement between two chains can cause 



their relative velocity to go to zero. For long enough 
chains, one would expect that entanglements would then 
dominate the collisions. However their effects usually 
only become dominant for chains larger than N « 500, 
as has been seen with the study of topological effects in 
flexible polymer rings [3] 

To understand the details of the dynamics of these col- 
lisions analytically is particularly difficult in light of the 
topological nature of the interactions and is beyond the 
scope of the present work. I would conjecture that the 
effects of entanglements will strongly suppress these for- 
ward scattering collisions for long enough chains. 

If the chains are of different sizes, then their behavior 
will depend on their relative sizes. From the above dis- 
cussion, the probability that a single monomer will collide 

— 1/5 

with a polymer of A'^i links is cx Therefore if a sec- 

ond polymer is of length N2 < Ni, the probability of one 
collision is approximately oc N2/nI^^ . In other words, 

1 /5 

the second molecule is unlikely to collide if iV2 ^ iVj^ . 
For chains of comparable sizes, the above results suggest 
that collisions will still be strongly inelastic and will also 
show behavior close to the thermal limit. The departure 
from this due to forward scattering will increase as the 
mass difference increases. 



IV. CONCLUSIONS 

In this work, collisions between polymers in a vacuum 
have been examined. This is quite an unusual situation 
because collisions typically occur between much smaller 
molecules, and consequently new features of this situa- 
tion are expected. The focus here is on chains where the 
initial relative center of mass velocity is so large that the 
initial internal thermal energy of the molecule can be ig- 
nored. This is an interesting limit to consider because 
for athermal chains, the magnitude of the initial velocity 
factors out. 

Typically, in a collision of such long molecules, they 



(a) 




(b) 




FIG. 6. (a) The probability density as a function of final speed v, and the amount of deflection z = cosO. (b) The probability 
density as a function of collision time t, and the amount of deflection z = cos 6. 



remain in contact for long enough to be close to equilib- 
rium so that they can be described by thermal averages. 
The final velocities and scattering angles, are close to 
such averages. However the main discrepancy with this 
description is due to a class of collisions that strongly for- 
ward scatter and remain in contact for a relatively short 
amount of time. This is evidence of strongly inelastic 
collisions of subsections of the chains. It is expected that 
for long enough chains, N > 500 entanglement effects 



will diminish the frequency of these kinds of collisions. 

In future work, it would be interesting to consider the 
case of chains with strong attractive forces as this may 
also be relevant to experiments and has interesting phys- 
ical consequences. For example, chains at temperatures 
that are initially below the coil-globule transition will 
tend to stick together for low enough relative velocities. 
However beyond a threshold velocity, they would be ex- 
pected to heat enough upon impact to become swollen 
and then separate from each other. 
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